Crossover from entropic to thermal dynamics in glassy models 



m 
o 
o 

(N 

S3 
S3 



S3 
S3 



I 

S3 
O 

o 



> 
00 

m 
in 
in 
o 
m 
O 

c3 



i 

S3 

o 
o 



X 
S3 



Eric M. Bertin 

CEA - Service de Physique de I'Etat Condense, 
Centre d'Etudes de Saclay, F-91191 Gif-sur-Yvette Cedex, France 
(Dated: February 2, 2008) 

We study the emergence of a crossover from entropically driven to thermally activated dynamics 
in different versions of the 'entropic' phase space model introduced by Barrat and Mezard, and 
previously considered in the zero temperature limit. We first focus on the low temperature (T <T g ) 
aging phase of the original model and show that a short time singularity appears in the correlation 
function for T g /2 < T < T g , leading to dynamical ultrametricity at T = T g . We then consider the 
finite size version of this model, showing that the long time dynamics is always thermally activated 
beyond a size dependent crossover time scale. We also generalize the model, introducing a threshold 
energy so as to mimic a phase space composed of saddles above the threshold, and minima below. 
In this case, the crossover time scale becomes much smaller than the equilibration time, and both 
kinds of aging dynamics are successively found, inducing a non trivial aging scaling which does not 
reduce to the usual t/t w (or even t/t^) one. 

PACS numbers: 75.10.Nr, 02.50.-r, 64.70.Pf 



A major part of our understanding of the glass transi- 
tion comes from Mode Coupling Theory (MCT) for struc- 
tural glasses 0,0 and from dynamical mean- field theory 
for spin glasses |3, H , which are the most sophisticated 
microscopic theories available in glassy physics. Apart 
from these theories, the energy landscapeparadigm has 
raised a lot of interest in recent years |a, la, [7|, |8| , due to 
the increasing power of computer simulations which now 
allow to implement the seminal ideas of Goldstein 0, 
specified later by Stillinger and Weber 0]. These stud- 
ies have led to associate the mode coupling transition 
temperature T c , where MCT locates an ideal glass tran- 
sition, with a threshold energy E t h of geometric nature, 
separating regions of phase space dominated by saddle- 
points for E > Eth from re gion s dominated by minima 
for E < E th [II 111 El El TOI In this framework, the 
ideal glass transition found in MCT and in mean-field 
theories is explained on the one hand by the sharp sep- 
aration between saddles and minima, and on the other 
hand by the divergence of the energy barriers surround- 
ing minima |l6j| . 

If one of these two conditions was not fulfilled, the 
ideal glass transition would be smeared out, which in- 
deed happens in real glass formers [2j. In particular, if 
energy barriers remain finite, one may expect not only 
the behaviour of the relaxation time with temperature 
to be smoothed, but also correlation functions to reflect 
this dynamical crossover between a time regime when 
saddles are mainly explored and a later regime in which 
the system visits minima in majority. The key point 
here is that visiting minima corresponds to thermally 
activated dynamics where the system has to overcome 
energy barriers, which thus depends crucially on temper- 
ature, whereas visiting saddles is associated to a rather 
temperature independent kind of dynamics, sometimes 
called 'entropic dynamics': the evolution of the system is 
mainly limited due to the lack of escape downwards di- 
rections, when visiting a saddle with only a few negative 



eigenvalues. Interestingly, this dynamical transition has 
already been observed in Lennard- Jones supercooled liq- 
uids [l7llla |. but a quantitative analysis of the associated 
crossover scales has not been achieved yet. 

In this paper, we study quantitatively this crossover 
and its dependence upon relevant parameters in simple 
stochastic models of glassy dynamics. We first consider 
the phase space model introduced by Barrat and Mezard 
[T^ | - hereafter denoted as BM model - which has known 
a recent surge of interest in the context of fluctuation dis- 
sipation relations pfl l2l| , and study the aging regime of 
the correlation function C (t w ,t w +t) in the whole glassy 
phase T < T g , generalizing previous results obtained in 
the limit T — > 0. We find a temperature independent long 
time behaviour (t 3> t w ), as well as a short time singu- 
larity (t <C t w ) in the temperature range T„ 



<T <T, 
for T = T, 



which leads to dynamical ultrametricity 
as in the trap model [23] ■ We then study the correlation 
function in the finite size BM model, both in equilibrium 
and in the aging regime, which allows to evidence clearly 
a dynamical crossover from entropic to activated dynam- 
ics. Finally, the BM model is generalized by introducing 
a threshold energy so as to mimic a phase space struc- 
ture composed of saddles and minima. We show that the 
change of aging dynamics (from entropic to thermally 
activated) in this model leads to a non trivial scaling be- 
haviour, and that the crossover time scale can be much 
smaller than the equilibration time. 



I. THE BARRAT-MEZARD MODEL 

One of the first attempts to describe aging dynamics 
with simple stochastic phase space models was the trap 
model proposed by Bouchaud (24|, with the aim to ev- 
idence a possible general scenario for glassiness, which 
emphasizes activation effects. A few years later, Barrat 
and Mezard |l9j | proposed a somewhat similar model, but 
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where thermal activation was replaced by an 'entropic 
activation', i.e. by the fact that at low temperature less 
and less lower energy states could be reached when the 
system is already in a low energy state. In a continuous 
energy formulation, corresponding to an infinite number 
of states, the BM model is defined as a Markovian pro- 
cess where a 'particle' (which usually corresponds to the 
representative point of a system in its phase space) is al- 
lowed to jump from one energy state to any other one. 
These energy states are distributed according to an a pri- 
ori distribution p(E) 1 and transition rates W(E'\E) from 
energy E to energy E' correspond to the Glauber choice: 



W(E'\E) 



T oP (E') 



l + e (E'-E)/T 



(1) 



where To is a microscopic frequency scale, that we shall 
set to unity in the following. The dynamics of the BM 
model is then described by the probability P(E, t) to have 
energy E at time t, which satisfies the following master 
equation: 

^ = y" dE> '{W(E\E')P(E> ,i) -W(E'\E)P(E,t)] 

• (2) 

where for simplicity - but without loss of generality since 
we are interested in the low energy dynamics - we have 
restricted p(E) to be non-zero only for E < 0. 



A. Zero temperature results 

Contrary to the trap model, and to any model driven 
by thermal activation, the temperature can be set to zero 
without freezing completely the dynamics, leading to an 
evolution somewhat similar to a steepest descent dynam- 
ics in a high dimensional continuous energy landscape. A 
noticeable property of this model in the zero temperature 
limit, which has been studied in details in |l9j . is that the 
dynamics, expressed in terms of sojourn times instead of 
energies, becomes completely independent of the func- 
tional form of the distribution p(E); in particular, p(E) 
can be bounded or not. Defining the average sojourn 
time t(E) at energy E by the relation: 



1 



r{E) 



dE'W(E'\E) 



(3) 



which reduces for T = to t(E)- 1 = /f^ dE' p(E'), one 
finds that the a priori distribution ipo( T ) derived from 
p(E) is given by ipo( T ) = t~ 2 9(t — 1), independently 
of p{E) . As a result the dynamical distribution p(r, t) , 
which has been computed in |l9j| , does not depend either 
on the shape of p(E). It was found that p(r,t) cannot 
reach a steady state, but exhibits for large t a scaling form 
p(r, t) = t _1 0o(r/i), typical of the aging behaviour, with 

<M«) s iven b y ; 



<po (u) 



-l/u 



(4) 



Defining the correlation function C(t w ,t w + t) as the 
probability not to move between t w and t w + t, which 
can be written: 

r° 

C(t w ,t w +t) = / dEPiE^e- 1 '^ (5) 

-' CO 

one finds for T — and large t w : 



t „ 



tin ~\~ t 



(6) 



Note that in the 'short time' regime 1 <C t <C t w , the 
correlation function C(t w , t w +t) behaves linearly, in the 
sense that 1 — C ~ t/t w ; this point will be important in 
the discussion of the finite temperature results. 



B. Aging at finite temperature 

If one chooses an exponential density of states p(E) = 
(3 g exp((3 g E) 9(—E), as is usually done in the trap model, 
then the Boltzmann equilibrium distribution becomes 
non normalizable for T < T g = l/{3 g . In this case, the 
dynamical distribution drifts towards deeper and deeper 
energies, leading to the aging phenomenon, which can be 
evidenced for instance by the behaviour of the correla- 
tion function. In this section, we present the new results 
we have obtained in the aging regime for the finite tem- 
perature case, namely for < T < T g . 

As was done for the trap model in [2q|. one can then 
look for a scaling solution for the dynamical energy dis- 
tribution. A natural scaling variable would be the ratio 
of the typical time r(E) the system spends in a state with 
energy E - as defined in Eq. J3J - to the time t elapsed 
after the quench. Computing t{E) for < T < T g and 
large \E\ yields: 



t(E) 



Sm7r M e -E/T g 
7T/1 



(7) 



where p = T/T g . We choose a scaling variable u propor- 
tional to r(E)/t, defining it through u^ 1 = te E / Tg . This 
leads to a distribution P(E, t) = j3 g u <fi(u), where <j>(u) is 
normalized by J Q °° du<ft(u) = 1. Once these substitutions 
are made in Eq. J2J), the following equation is obtained: 



u 2 ct>'{u) + (u-J) (j>{u) 



dv 



1 + (v/u) 1 /^ 



(8) 



where J = np,/ sin 7174. The asymptotic behaviour of the 
distribution P(E, t) can be extracted rather easily from 
this equation. Let us analyse for instance the large u 
behaviour of <j>{u) - or large \E\ behaviour of P(E, t). In 
this case, the integral in the r.h.s of Eq. (JSJ reduces to: 



00 Hv) , 

dv . — , 

o 1 + (v/u) 1 /^ J 



dv <j>(v) = 1 u > 1 (9) 
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FIG. 1: Plot of the aging correlation function C(t w ,t w +t) of 
the BM model, for t w = 10 4 , 10 s , 10 6 and 10 7 (from left to 
right), and temperature /i = 0.7. Inset: correlation plotted 
as a function of the rescaled time t/t w , showing a very good 
collapse of the data, apart from small deviations for t <C t w 
due to finite time effects. 



so that Eq. JSJ reads, neglecting also J with respect to 
u and calling </>/ (u) the asymptotic large u expression of 
0(«): 



u 2 (j>[(u) + u (f>i(u) = 



The general solution of this equation reads: 

A 1 
- + — 

11 11-^ 



4>i{u) 



(10) 



(ii) 



The constraint that 4>(u) must be normalizable requires 



that A — 0, so that 4>i(u) 



Interestingly, we find 



here the exact asymptotic expression for (/)(u), including 
the prefactor - equal to 1 - although we performed only 
an asymptotic analysis on a linear equation. Note also 
that taking into account the leading correction to the 
approximation made in Eq. @ leads to 4>i(u) — uT 2 ~ 
it~ 3 . Neglecting these corrections, one thus finds that 
P(E, t) is equal to f3 g t e E l T<1 - i.e. proportional to p(E) 
- for el £ l/ T « » t. 

In the opposite case »<1, the equation governing the 
asymptotic expression 4> s (u) of 4>(u) reads: 



u 2 4>' s (u) - J(j) s {u) = 



dz 



4>s{uz) 



Searching for a power law solution (t> s {u) 
above equation is satisfied only if: 



a u 



J 



dz 



1 + z 1 /^ 



(12) 



the 



(13) 



since the first term u 2 (j)' s (u) can be neglected in this case. 
Replacing by the respective values of both sides, one has: 



The convergence of the integral appearing in Eq. i|13|) 
requires that < (1 + a)p < 1. In this range, Eq. (|14|) 
admits two solutions: a = and a = l//i — 2. Since 
the zero temperature limit of the model is regular, the 
finite temperature solution must converge towards the 
zero temperature distribution given by Eq. (@J), which 
decays faster than any power law for u — > 0, so that one 
expects a — > oo for /i — > 0. The correct solution for (j){u) 
is thus: 



0(u) 



l/p-2 



U < 1 



So the energy distribution behaves at small \E\ 

e \E\/T g ^ t _ ag . 

P(£,i)~t 1 - 1 //' e (l-l//0iS/T. 



i.e 



(15) 
. for 



t 1 - 1 ^p(E)e- E / T 

(16) 

which means that small energies are almost equilibrated. 

Turning to the correlation function C (t w ,t w + t) de- 
fined above, Eq. (0) can be rewritten in terms of the 
scaling function 4>(u), using also Eq. J7J): 



C(t v 



t) 



dvcf>{v)e- Jt/vK 



(17) 



which shows that the correlation function exhibits a sim- 
ple aging form: 



C(t W: t w -\- £) — C a g(t j ^t w ^ 



(18) 



Using the above asymptotic results for (f>(u), one can find 
the short time and long time behaviour of the correlation 
function. In particular, for times t^>t w , C(t w , t w + 1) is 
given by: 



sin7r/i t w 

7T/i t 



t > u 



(19) 



for any /i < 1, so that the 'tail' of the correlation func- 
tion appears to be a temperature independent power law, 
apart from the prefactor. This point will be important 
for later discussions. 

The short time behaviour is computed from: 



1 C[t w , t w + t) 



dv(l)(v){l-e- Jt/vt ™) (20) 



For t <^ t w , one can try to linearize the exponential in the 
last factor. Using Eq. l[T5)l. one sees that the resulting 
integral converges at its lower bound only if /i < |; in 
this case, the rescaled correlation function C ag (u) is then 
regular when u — > (i.e. 1 — C ag (it) ~ u), just as in the 
zero temperature case. On the contrary, for | < /j, < 1, 
a singularity appears, and one then finds: 



1-C a 



(x-n)/n 



t < U 



(21) 



This property can be interpreted as follows: the average 
energy (E') E reached in a transition from E to E' reads: 



7T/Z 



7T// 



sin7r^ sin(l + a)7r/i 



(14) 



(E') 



dE' E'W{E'\E) = E 



TTfi 



tan 7171 



T g (22) 
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FIG. 2: Short time behaviour of the correlation function; 
1 — C(t w , tw+t) is plotted as a function oft/t w for several tem- 
peratures fx = 0, 0.2, 0.4, 0.6, 0.7 and 0.8 (from right to left), 
showing the onset of a short time singularity for i < fx < 1 
with an exponent (1 — fl)/fi (dashed lines). For fx < 1 — C 
is linear for t <C t w . Inset: long time tail of the correlation 
for /j, = 0.2, 0.4, 0.6 and 0.8 (from right to left), showing a 
power law decay C ~ (t/tm) -1 . 



For /j < i (E') E < E so that the dynamics resembles the 
zero temperature one: at each step, the energy is lowered 
on average. On the contrary, for /i > |, (E') E > E and 
the energy is raised on average at each transition. So 
the slow drift of the average energy E(t) as a function of 
time results in this case from a kinetic effect: the particle 
stays a longer time on lower energy states. Actually, the 
particle typically jumps a lot of times among high energy 
states before reaching a low energy one, then jumps back 
to high energies and so on. This scenario is reminiscent of 
what happens is the trap model, for which (E') E = T g in- 
dependently of the starting energy E, and is presumably 
responsible for the onset of the short time singularity. 
Note however that in the BM model, both entropic and 
thermal dynamics come into play for 1/2 < fx < 1, the 
latter being rather irrelevant for ll < 1/2, but thermal ac- 
tivation remains in any case an auxiliary mechanism, the 
main one being entropic. Indeed, this can be seen from 
the typical sojourn time t(E) - see Eq. (jTJ) - which is al- 
ways proportional to e~ E ^ Ta , and not e~ E l T as would be 
the case for a purely thermal dynamics. So the present 
scenario is a bit different from the usual landscape pic- 
ture, since the entropic channel exists here even at very 
low energies. We shall see in the next sections how to 
generalize the model in order to account for a complete 
vanishing of downwards escape directions. 

We now present numerical results in order to confirm 
the above analysis. For computational reasons, simula- 
tions were performed with a slightly different version of 



the model, in which the transition rates are given by: 

W(E'\E) = T p(E') E'<E (23) 

W(E'\E) = T Q p{E')e- {E '- E ^ T E' > E (24) 

These transition rates W(E'\E) behave qualitatively as 
the usual Glauber rates, so that one expects this differ- 
ence not to affect the scaling and power law behaviours, 
but only the numerical prcfactors. One needs first to 
check the aging form - Eq. I|18fl _ of the correlation func- 
tion. This is done on Fig. ^ which displays the correla- 
tion function C(t w , t w + t) for /i = 0.7 and several waiting 
times, t w = 10 4 , 10 5 , 10 6 and 10 7 . The rescaling, shown 
in inset, is very good, with only small deviations in the 
short time regime t <C t w where finite time effects be- 
come stronger. We also test the predicted exponents for 
the short time and late time behaviour of the correlation, 
as shown on Fig. EI One clearly sees the linear behaviour 
of 1 — C for /i < i, as well as the predicted singularity 
with exponent (1 — fx)/ fi for i < fx < 1. It is however 
difficult to find a genuine straight line for this short time 
singularity since t must satisfy 1 <C t <C t w , and such a 
separation of scales is hard to obtain numerically (note 
that t w is already as large as 10 7 ). The late time regime is 
easier to evidence quantitatively, since the only condition 
is i 3> t w ; the corresponding power law behaviour pro- 
portional to (t/tyj) -1 , with a temperature independent 
exponent, is clearly evidenced on the inset of Fig.0 As 
explained above, we cannot test quantitatively the pref- 
actor predicted in Eq. I|19|) since the numerics deals with 
a slightly different model, but we can see that the prefac- 
tor found in simulations decreases with temperature, in 
qualitative agreement with the prediction. 



C. Dynamical ultrametricity for T = T g 

The study of the behaviour of the model right at the 
glass transition temperature T — T g appears to be easier 
to solve, although the scaling form is a bit more subtle. In 
order to simplify the notations, we choose in this section 
T g as the energy (and temperature) unit. It can be shown 
that the correct scaling variable u is given in this case by 
u^ 1 = e E tint, and that the associated scaling function 
<p(u) must be related to P(E, t) through: 



P{E,t)~— u<p(u) 



(25) 



The normalization condition f_ dE P(E,t) translates 
into: 



- — / du(p(u) 



lnt 



1 



t — ► oo 



tint)" 



(26) 



which requires that tp(u) ~ u 1 for u — > 0. The master 
equation then becomes in the limit of large t: 



u if'(u) + (u — 1) <p(ix) 



1 

u 



(27) 



5 



The condition that <p(u) must be normalizable in the 
sense of Eq. (|26[1 selects a unique solution, which reads: 



1 



p( u ) = -(l_ e -V») 



(28) 



so that P(E, t) is given by 



1 



P(E, t) = — [l- exp(-e E t In t)] (29) 

From this expression, one can compute the correlation 
C(t w ,t w + t) at the critical temperature: 

C{t w ,t w +t)= dEP{E,t w )e- t/T{E) (30) 



with t(E) given by 
1 



t(E) 



e E ' dE> 



1 + e 



E'-E 



e E hx{l + e- E ) (31) 



which reduces to \E\ e E for large \E\ (E < 0). So the 
correlation reads, for large t w and i: 



C(t w , t w +t) = 



1 



lnt 



o 



dEexp{Ete E ) (32) 



X 



[l - exp(-e E t w lnt w )] 



From this expression, one can deduce that the rele- 
vant scaling variable is no longer t/t w , but rather to — 
hit/ hit w , in the sense that when taking the infinite t w 
limit keeping u> fixed, C(t w , t w +t) converges to a function 
C(ui) given by: 

C(u)) = l-u (0<u<l), C(u)=0 (w>l) (33) 

This is precisely the asymptotic form found in the trap 
model at the glass transition temperature, which was 
shown to satisfy dynamical ultrametricity [23, u3 • So we 
see that dynamical ultrametricity is indeed independent 
of the type of dynamics considered - entropic or thermal 
activation - but is rather governed by the exact balance 
of the Boltzmann weight and of the exponential density 
of states which is valid only right at the glass transition 
temperature. Interestingly, the differential equation - 
Eq. H27|) _ satisfied by the scaling function (p(u) appears 
to be the same in the BM model as in the trap model 
- although another approach was used in |2^| - lead- 
ing to the same scaling functions </?(.), even though they 
are derived from two different master equations. Note 
however that the scaling variables are slightly different, 
u = {t\nte E y 1 in the BM model, and v = {te E )~ 1 in 
the trap model. 



II. FINITE SIZE EQUILIBRIUM 

In realistic models of glasses, the structure of the global 
phase space is expected to be highly non trivial, since it 



has to encode the finite dimensionality of the physical 
space. So it appears tempting, at least as a first approx- 
imation, to think of a macroscopic system as a collection 
of small independent subsystems, the size of which would 
be of the order of the coherence length |26l l27| - see also 
28] for an experimental investigation of this point. This 
scenario has indeed been supported by recent simulations 
|29| on Lennard- Jones supercooled liquids, showing that 
a system with 130 particles has to a good approximation 
the same dynamical behaviour as two non interacting 
systems composed of 65 particles each [3^ . 

In this context, finite size versions of the usual - trap 
and BM - phase space models appear to be particularly 
relevant, and the simplest case to study is the equilibrium 
one. In order to make contact with more realistic models, 
let us emphasize that 'states' in the present case should 
be associated with singular points in phase space, namely 
saddles or minima, so that one expects their number M 
to grow exponentially with the number N of physical 
degrees of freedom: 



M ~ M e 



(34) 



where a is a numerical constant. A lot of studies have 
investigated the value of a as a function of temperature 
when considering only minima and have found a -C 1 
El • Note however that, bearing in mind the above 
picture, N should be itself a function of temperature since 
it depends on the coherence length £(T) through N ~ 

Considering M energy states randomly chosen from 
an exponential distribution p(E), the equilibrium dis- 
tribution concentrates onto the lowest energy states, so 
that the correlation function decays on a time scale 
corresponding to the typical largest sojourn time r max : 
C eq (t,M) — C(t/r max ). The value of r max can be esti- 
mated using a scaling argument: in a finite size model, 
one can imagine that all the sojourn times are chosen 
at random according to the probability distribution ip(r) 
associated to the infinite size model. In the exponential 
trap model, ^ tm (t) = /j/t 1+m and this picture is indeed 
exact, since the relation r = e E l T holds locally, at the 
trap level: one has then r max ~ M 1 ^. On the contrary, 
in the BM model, the relation given in Eq. @ is valid 
only in the infinite size limit, and should be only an ap- 
proximation in the finite M case. Still, one can hope that 
this approximation is able to give the correct scaling of 
T max with M: ip(T) ~ 1/r 2 then yields r max - M. 

An interesting question naturally arises about the fi- 
nite size equilibrium dynamics in the low temperature 
phase. Given that for T < T g the continuous energy 
equilibrium distribution P eq (E) oc p{E) e~ E / T becomes 
non normalizable, a cut-off scale necessarily appears in 
order to regularize the distribution P eq {E). One way to 
introduce such a cutoff is precisely to consider a finite 
number of states. Otherwise, the spontaneous cutoff is 
generated by the dynamics itself, since it is given by the 
time t w elapsed after the quench. So one can wonder 
whether or not the dynamics is qualitatively the same 



0.8 



0.4 
0.2 



1 1 












1 

sr °- 8 


5^ 1 


— W- 






— \ \ 


tt tr 0.6 










\ 


\ \ °- 2 
\ \ \ 




1 . 








10 2 


10° 


10 2 








t/M 




1 , 











^ I — i i ~- =r--4-^= ^-f-^- 1 

10° io 2 io 4 io 6 io 8 io 10 



I 

FIG. 3: Plot of the finite size equilibrium correlation function 
C eq (t,M) of the BM model, for M = 10 3 , 10 4 , 10 5 and 10 6 , 
and temperature fi — 0.7. Inset: correlation plotted as a 
function of the 1 time t/M, showing a good collapse of the 
data, except for short times where the M = 10 3 curve does 
not rescale completely. 

whatever the nature of the cut-off; we shall see that the 
answer depends on the model considered. 

In the exponential trap model in the low temperature 
phase fi < 1, with a finite number M of traps, the cor- 
relation C™(t, M) is thus a function C T ^(t/M x /< Jl ) of the 

rescaled time t/M 1 ^, which can be computed in the 
asymptotic regimes t <§; M 1 ^ and t M 1 ^, leading 
to the interesting result that the behaviour is similar to 
that found in the aging case, replacing the cut-off t w by 
M 1 /^. However, the prefactors are different from that 
found in the aging case [32|. Namely, one finds: 

l-C™(t,M) ~ n s (j^j " t«M 1 /M (35 ) 
where the prefactors n a and ki are given by: 

Ka = -^-r( / i)r(i/ / i)r(i-/x)- 1 //' (37) 
l-fj, 

Kl = %/2 1_2M Alv / ^rQ + ^ (38) 

This shows that in the trap model, the asymptotic be- 
haviour of the rescaled correlation function for ji < 1 is 
the same for t much smaller or much greater that the cut- 
off time scale, whatever the origin of this scale, namely 
static (number of sites) or dynamic (waiting time). It 
would be interesting to see if this property also holds for 
the BM model. 

To this aim, we first check numerically the proposed 
scaling in t/M for the BM model, using in this case 
the natural version of the model, with Glauber rates. 
Fig. |3 displays C eq (t, M) for several values of M, namely 




t/M 



FIG. 4: Asymptotic large time behaviour (t 3> M) of the cor- 
relation Ceq(t,M) with M = 10 3 for temperatures /i = 0.2, 
0.3, 0.5 and 0.7 (from top to bottom), showing a tempera- 
ture dependent power law tail C sq (t, M) ~ (t/M)~~^ (dashed 
lines), typical of activated dynamics. Inset: short time be- 
haviour (t < M) of the correlation, with M = 10 6 , /U = 0.3, 
0.4, 0.6, 0.7 and 0.8, showing the onset of the same short 
time singularity as in the aging regime (see Fig. dotted 
line: slope — 1). 

M = 10 3 , 10 4 , 10 5 and 10 6 ; the inset shows that the data 
collapse well when plotted as a function of t/M, with 
again small deviations in the short time regime which re- 
quires larger system sizes to reach the asymptotic scaling 
form. 

Turning to the asymptotic behaviour of the correla- 
tion, one sees on Fig. 0] that the late time behaviour 
of C eq (t,M) is no longer temperature independent, as 
in the aging case, but depends strongly on temperature, 
and the slope is precisely equal, to numerical precision, 
to — ji. This is a strong evidence for a thermally driven 
dynamics at large times, which corresponds to the fact 
that the particle can only escape from the deepest traps 
through thermal activation. On the contrary, the short 
time singularity is not modified with respect to the ag- 
ing case, as can be seen in the inset of Fig. This was 
expected as the short time dynamics is dominated by 'in- 
termediate' energy states - i.e. the highest states among 
the low energy ones onto which the equilibrium measure 
concentrates - from which one can escape downwards. 

As a result, a crossover appears in the equilibrium fi- 
nite size correlation between an entropic behaviour at 
short times and an activated one at long times. A natural 
question at this stage would be to see if this crossover is 
a peculiarity of the equilibrium case, or if it can be found 
also in dynamical regimes. 

III. CROSSOVER IN THE AGING REGIME 

Thermally activated phenomena are known to play an 
important role in glassy physics, and numerous experi- 
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FIG. 5: Correlation C(t w ,t w +t, M) in the out of equilibrium 
regime of the finite size model, for M = 10 3 , t w = 10 2 and 
u. — 0.2, 0.3, 0.5 and 0.7. One sees a crossover from a rather 
temperature independent (entropic) regime to a temperature 
dependent (activated) one with exponent —u, (dashed lines). 
Inset: Effect of a Unite connectivity Z, for M = 10 3 , t w = 10 2 
and [i, = 0.5; dotted line: Z = 100, dashed line: Z = 300 
and full line: fully connected case. Reducing the connectivity 
decreases the crossover time scale. 



mental studies have focused on the effect of temperature 
on glasses (see e.g. 33] for structural glasses and [3^ for 
spin glasses). Theoretical works also support this view, 
since the ideal glass transition predicted by MCT has 
been shown to be smeared out in real glasses by acti- 
vated processes 0, 0. However, in structural glasses, 
these activated processes are usually assumed to be rel- 
evant only on time scales close to the equilibration time 
t(T) - which becomes very large below the glass tran- 
sition temperature - and not in the aging regime which 
corresponds to shorter time scales. Indeed, recent simula- 
tions on Lennard- Jones liquids seem to validate a mean- 
field-like scenario with a well-defined effective tempera- 
ture [3{| , at odds with the behaviour of thermally driven 
models like the trap model j3|| - except in some very 
special cases pflj . 

One might thus think that the crossover from entropic 
to thermal dynamics observed in numerical simulations 
0>H] is om y due to the small size of the systems consid- 
ered, and that the associated time scale should be very 
large for macroscopic systems, as shown in the mathe- 
matical analysis of the dynamics of the random energy 
model [221 E3- Still, the arguments presented above 
suggest that for a finite dimensional system, this time 
scale should be governed by the coherence length rather 
than the global size of the system. One can also wonder 
whether the crossover time scale could be possibly much 
smaller than the equilibration time under some condi- 
tions. 

In this section, we study the onset of a dynamical 
crossover from entropic to activated behaviour in the ag- 
ing regime of the correlation function. To this aim, we 



consider the finite size BM model as well as a gener- 
alization including a threshold energy below which all 
states are isolated minima, and show that the dynamical 
crossover time is equal to the equilibration time only for 
the fully connected model, and becomes much smaller as 
soon as a non trivial phase space structure is taken into 
account. 



A. Aging in a finite size system 

A natural generalization of the analysis performed in 
the previous section is to study the aging regime 1 <C 
t w <C M of the finite size BM model, with the aim to 
evidence a crossover in the correlation C(t w ,t w + t,M) 
from the usual BM aging form for t <C M to an activated 
behaviour for t 3> M, So the correlation function is 
now characterized by a couple of time scales, t w and M, 
instead of only one of these. 

Numerical results for the correlation function 
C(t w ,t w + t,M) are shown on Fig. for M = 10 3 , 
t w = 10 2 and various temperatures (/1 = 0.2, 0.3, 0.5 
and 0.7). Apart from a global shift, the 'early parts' 
of the curves (i.e. times t such that t w < t < M) are 
almost temperature independent, which is characteristic 
of an entropic kind of dynamics, as we have seen in the 
previous sections. On the contrary, the late [t > M) 
power law decay, with an exponent —fx (dashed lines) is 
the fingerprint of thermally activated dynamics. 

As we have only considered up to now fully connected 
models, one can wonder whether a finite connectivity 
Z <C M would strongly or not modify the previous re- 
sults. Although a detailed study of the BM model with fi- 
nite connectivity is beyond the scope of the present work, 
we comment here on the results of preliminary simula- 
tions on random graphs with finite average connectivity. 
In particular, decreasing the connectivity does not af- 
fect much the overall decay time of the correlation func- 
tion, in equilibrium as well as in the aging regime. Still, 
out-of-equilibrium numerical data presented on the inset 
of Fig. for M = 10 3 and different values of Z, indi- 
cate that reducing the connectivity tends to decrease the 
crossover time scale beyond which thermal activation be- 
comes dominant. This could have been expected, since a 
finite connectivity generates a lot of local minima, from 
which one can escape only through thermal activation. 
So we can guess that in a realistic glass former, even 
though the 'elementary correlated cell' can contain a lot 
of degrees of freedom, the effect of connectivity can fur- 
ther reduce the crossover time scale, suggesting that ac- 
tivated phenomena should be relevant on accessible time 
scales. 

Although a quantitative study of the finite size BM 
model in the out-of-equilibrium regime would be of great 
interest, it would require a larger scale separation than 
the one we obtained, which is hard to reach for com- 
putational reasons. We thus propose another version of 
the BM model, in which one can safely take the limit 




t 



FIG. 6: Plot of the correlation function C(t w , t w + t) of the 
ThBM, in the aging case (/i — 0.7) and in the regime t^>t w 
(long time) and t w <C r t h- Both r t h and t w are varied: Tth = 
10 5 (t w = 10 2 , 10 3 j, 10 6 and 10 7 (t w = 10 2 , 10 3 , 10 4 j. A 
crossover appears between an entropic regime, with exponent 
— 1, and an activated one, with exponent —fi. 

M — > oo without eliminating the dynamical crossover we 
are presently interested in. This can be done by taking 
into account a non trivial structure of phase space. 

B. Effect of a threshold energy 

Coming back to the paradigm of the energy landscape 
with a sharply defined separation in energy between sad- 
dles and minima [Til [13 . PHI , Il4l TTEj . one could try to 
mimic this phase space structure by introducing by hand 
a threshold energy E t h acting as a kinetic constraint, 
such that direct transitions between states with energies 
E and E' are forbidden if E and E' are both below E t h, 
and otherwise unchanged with respect to the usual BM 
model. In other words, states with energy E > E t h can 
be considered as saddles, from which all other states are 
accessible, whereas states with energy E < E t h repre- 
sent isolated minima from which the system can escape 
only by jumping to a saddle, before eventually reaching 
another (possibly deeper) minimum. Such a definition 
could be applied to a finite size model, and the equili- 
bration time would be in this case a function of M. In 
Lennard-Jones liquids, this threshold E t h has been found 
to be proportional to the number N of degrees of freedom 
but, as argued above, N should be actually bounded by 
£(T) d , so that Eth is expected to remain finite even for 
macroscopic systems. Conversely, M should also remain 
finite, but assuming that the smallest energy E m i n found 
in the system lies much below E t h , one can still take the 
limit M — ► oo keeping Eth fixed, so as to describe the dy- 
namics on time scales much smaller than the finite size 
equilibration time. 

Fig. H3 displays the results of numerical simulations of 
the correlation function in this generalized model, de- 




FIG. 7: Same curves as on Fig. [S| rescaled in their late 
part in order to evidence the crossover scaling, by plotting 
(Tth /t w ) C(t w , t w + t) versus t/r t h (see text). Dashed lines 
indicate slopes —1 and —fi respectively. Inset: rescaled cor- 
relation function for [i = 0.2, 0.3, 0.5 and 0.7 (from top 
to bottom), t w — 10 3 and Tth = 10 7 , showing the strong 
temperature dependence for t 3> Tth, whereas the slope for 
t w <C t -C Tth is independent of fj,. 

noted hereafter as ThBM ('Threshold Barrat-Mezard') 
model, for several values of t w and Eth, focusing on the 
long time behaviour (t t w ). The simulations show the 
onset of a crossover between two different power law be- 
haviours, the first one with an exponent —1 correspond- 
ing to the usual long time behaviour of the BM model 
(entropic dynamics) , and the other one with an exponent 
—\i typical of the long time tail of the correlation in the 
trap model (thermally activated dynamics). 

Let us now present a scaling argument in order to gain 
a better understanding of this crossover. We first focus 
on the most interesting case for which the typical ener- 
gies reached at t w , of the order of — T g lnt w , remain far 
enough above E t h- 

We assume, and we have checked numerically, that the 
tail (E much below — T g hit w ) of the dynamical energy 
distribution P(E,t w ) keeps the usual form P(E,t w ) oc 
e E / T a found both in the BM and in the trap model, for 
E > E t h as well as for E < E t h- The natural time versus 
energy relation in the BM model allows to define a time 
scale Tth associated to E t h through Tth = ex P{\Eth\/T g ) - 
we neglect here the prefactor appearing in Eq. Q. One 
thus expects the following asymptotic behaviour for the 
dynamical distribution p(r, t w ) of sojourn times: 

p{T,t w ) ~ i (M T<T th (39) 

tw \ T ) 

A ft \ 1+AI 
P(r,t w ) ~ r f-^J r»r t/l (40) 

these expressions corresponding to the large r tail of the 
dynamical distribution in the usual BM model and in 
the trap model respectively. The rescaling factor A is 
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FIG. 8: Short time behaviour of C(t w ,t w + t) for fj, = 0.3 
and t w » Tth, plotted as (t,„/r th ) 1_M [1 - C(t w ,t w + t)] versus 
t/rth [see Eq. J44D / so as to evidence the crossover from a 
linear (entropic) behaviour for t <C Tth to an activated one, 
with exponent 1 — fx, for Tth t <C t w ; r t h = 10 3 (t w = 10 5 , 
10 6 ; and 10 4 (t w = 10 6 , 10 7 ). Inset: same plot with fi = 0.6, 
such that the entropic regime 1< f < Tth corresponds to a 
power law with exponent (1 — 1^)/^; T t h = 5.10 3 , 10 4 , 2.10 4 
andt w = 10 s , 10 6 , 10 7 . 

determined by matching the two asymptotic expressions 
for r = Tth, which yields A = (£m/T t /i) 1_M . Computing 
the correlation for t <C T t h, one recovers as previously 
C(t w ,t w + i) oj t w /t, whereas for t > Tth, one now has: 

C(t w ,t w +t)~hL (IfY (41) 

Tth V t / 

Summarizing the above results, one can introduce a scal- 
ing function /(.) such that C(t Wl t w +t) reads: 

C(t w ,t w +t) = ^f(— ) (42) 

Tth \Tth J 

with the following asymptotic behaviour for f(z): 

/(*)-- (z>l) (43) 

z z^ 

In order to test this prediction, we have plotted on 
Fig. [7| the previous curves using the rescaled variables 
{Tth/tw) C(tw,t w + t) and t/Tth- The resulting collapse 
of the late part of the curves (i.e. for t ^> t w ) is very 
good, confirming the validity of the above scaling anal- 
ysis. Note that the whole curves cannot be rescaled in 
a simple way since the early part of the correlation still 
scales in t/t w as usual (rescaling not shown), so that two 
different scaling regimes t ~ t w and t ~ T t h have to be 
considered separately. Fig. [7| shows that in the limit of 
a large separation of scales between t w and T t h, the plot 
would reduce to two straight lines with different slopes 
— 1 and —fi. The temperature dependence is clearly evi- 
denced in the inset of Fig.[7| showing as expected that the 
slope is independent of fi for t w <^ t <^ T t h (and equal to 



— 1), whereas it is strongly temperature dependent (and 
equal to — /i) for t T t h- 

In the opposite regime t w ^> T t h, the system is trapped 
into deep minima so that the behaviour becomes similar 
to that of the trap model, up to an energy shift. However, 
the particle necessarily visits states with energy E > E t h 
in order to escape from minima, and then relaxes towards 
another minima. So a crossover between entropic and 
activated dynamics is also expected in this case for t ~ 
T t h, which now falls in the short time regime (1 <C t <C 
t w ). In order to analyse this crossover, one can derive 
a scaling relation along the same lines as for Eq. 142fl . 
which now reads: 

i-c(t w ,t w + t) = h^\ " g(^A (44) 

where the scaling function g(z) has the following asymp- 
totic behaviour: 

g{z)~zi (z<l), g{z)^ z 1 -^ (z > 1) (45) 

with 7 = min[l, (1 — A*)/mL as f° r the usual Barrat- 
Mezard model. Fig. [S] displays the numerical results for 
different t w and Tth, with ll = 0.3 and 0.6, and shows 
a good collapse of the data, which clearly supports the 
existence of a crossover for t ~ Tth- 

In summary, a non fully connected phase space struc- 
ture generates a large number of local minima, so that the 
dynamical crossover time scale associated to the end of 
the entropic regime is strongly reduced with respect to 
the fully connected case, suggesting that this crossover 
may be relevant in realistic glasses on experimental time 
scales, even below the glass transition temperature. 

IV. CONCLUSION 

In this paper, we have shown that a crossover from 
entropic to activated behaviour arises naturally in glassy 
models once a finite coherence length is taken into ac- 
count, allowing to conceptually decompose a large sys- 
tem into a collection of independent small subsystems 
[26l 1271 . In this framework, one is then naturally led to 
consider finite size phase space models to describe the 
local dynamics within an 'elementary correlated cell'. 

Interestingly, calculations in the trap model show that 
the finite size M replaces in a rough sense the time cut-off 
t w appearing in the aging regime, without affecting the 
asymptotic power law short and late time behaviour of 
the correlation. On the contrary, in the BM model, the fi- 
nite size equilibrium behaviour differs significantly from 
the aging one: the short time singularity remains the 
same, but the long time one is no longer governed by an 
entropic mechanism, but corresponds to a thermally ac- 
tivated dynamics since escaping from the deepest states 
is only possible through jumps to higher states, if the 
number of states is finite. Generalizing this analysis to 
the aging regime of finite size systems, where both t w 
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and M come into play (t w <C M) , the previous crossover 
still exists for t ~ M as in equilibrium, and can even be 
found earlier if a finite connectivity is taken into account. 
Although such a crossover has already been observed in 
numerical simulations |171 Il8| , a quantitative analysis of 
the associated time scales would be very valuable, in or- 
der to see how it varies with temperature and with the 
total size of the system. In particular, if a decompo- 
sition into small correlated cells holds, one expects the 
crossover time scale to saturate to a finite value for large 
system sizes. 

Following the phase space decomposition into saddles 
and minima, we have then studied a generalization of 
the BM model in which the crossover time scale is not 
induced by the finite size of the system, but instead by a 
threshold energy E t h which separates saddles for E > E t h 
from minima for E < E t h- Assuming that the equilibra- 
tion time is large enough, we find that the correlation 



function keeps aging even beyond the crossover time scale 
Tth, but with a behaviour characteristic of activation dy- 
namics. Interestingly, the transition from one type of 
aging dynamics to the other leads to an unusual scaling, 
due to the fact that two time scales t w and Tth are now 
involved instead of one. Moreover, the presence of a lot 
of local minima in phase space leads to a crossover time 
scale much smaller than the equilibration time, contrary 
to what happens in the fully connected case. 
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